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Introduction 

T HE inclusion of CFD-based analyses into dis- 
ciplines such as aeroelasticity, aeroservoelas- 
ticity, and optimization is currently not routine 
due to the high computational costs of the CFD 
analyses. One solution to this problem is the 
development of reduced-order models (ROMs). 
Reduced-order models capture the essence of the 
dynamical system under investigation, resulting in 
a significantly less complex model. The reduced 
complexity of the ROM translates into improved 
computational efficiency. The ROM can then be 
used for subsequent analyses at reduced computa- 
tional costs. It should be stated that the ROMs 
discussed within this context are a generalization 
of the traditional ROMs that involve the reduc- 
tion of the dimensions of a given matrix. So, for 
example, the lift response of a three-dimensional 
aeroelastic CFD system undergoing plunging mo- 
tions can be characterized, as will be shown, using 
a single impulse response function that relates lift 
due to plunge. In this case, the single impulse re- 
sponse function is the ROM. Therefore, a ROM 
can also be a functional condensation of the origi- 
nal system in addition to the traditional interpre- 
tations with respect to matrices. 

Early mathematical models of unsteady aerody- 
namic response capitalized on the efficiency and 
power of superposition of scaled and time shifted 
fundamental responses, or convolution. Classical 
models of two-dimensional airfoils in incompress- 
ible flow 1 include Wagner’s function (response to 
a unit step variation in angle of attack), Kuss- 
ner’s function (response to a sharp-edged gust), 
Theodorsen’s function (frequency response to sinu- 
soidal pitching motion), and Sear's function (fre- 
quency response to a sinusoidal gust). As geomet- 
ric complexity increased from airfoils to wings to 
complete configurations, the analytical derivation 
of response functions was no longer practical and 
the numerical computation of linear unsteady aero- 
dynamic responses in the frequency domain be- 
came the method of choice . 2 * And, when geometry - 
and/or flow-induced nonlinearities aerodynamic ef- 
fects became significant, the nonlinear equations 
were computed via time integration. The direct 
time-integration approach for solving aeroelastic 
problems via the coupling of the nonlinear aerody- 
namic equations with the linear structural equa- 
tions has yielded a very powerful simulation ca- 
pability with two primary challenges. The first 
challenge is the computational cost of this simu- 
lation, which increases with the fidelity of the non- 
linear aerodynamic equations to be solved. The 
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second challenge is that simulations cannot be used 
effectively for preliminary aeroservoelastic design. 
Design by simulation inevitably becomes design 
by trial-and-error, an impractical approach. The 
development of ROMs is targeted precisely at ad- 
dressing these two challenges. While a ROM pro- 
vides increased computational efficiency over the 
original more complex system, it is the mathe- 
matical model that is extracted from the original 
system that enables t he interconnection with other 
disciplines. This simplified reduced-order mathe- 
matical model transforms a highly complex ’’black 
box into a system with distinct physical and 
mathematical properties suitable for design analy- 
ses. Although the first challenge can be improved 
via parallel processing and advanced algorithms, 
the direct time-domain approach does not address 
the second challenge. 

Attempts to address t he problem of high compu- 
tational cost include the development of transonic 
indicia! responses . 3- 5 Transonic indicia! (step) re- 
sponses are responses due to a step excitation of 
a particular input , such as angle of attack, about 
a transonic (or nonlinear) steady state condition. 
More recently, Marzocca et al 6 have analytically 
derived indicia! functions for three-dimensional 
wings in compressible flow. Neural networks have 
also been used to develop nonlinear models of un- 
steady aerodynamics 7 and nonlinear models of ma- 
neuvers (using an experimental database ). 5 Neural 
networks and Volt, err a series have some similari- 
ties. since each involves the characterization of a 
system via an input-output, mapping . 7,8 In par- 
ticular, there is a direct relationship between the 
weights of a neural network and the kernels of a 
Vo] terra, series representation for a particular sys- 
tem . 9 

A major difference between Volterra series and 
neural networks, however, is in the training effort. 
Neural networks can require a substantial train- 
ing effort 0 while Volterra series require neither a 
training period nor curve fitting for model con- 
struction. Also, Volterra kernels provide a direct 
means for physical interpretation of a system’s re- 
sponse characteristics in the time and frequency 
domains. However, potential disadvantages of the 
Volterra theory method include input amplitude 
limitations related to convergence issues and the 
need for higher order kernels . 10 

Another approach for reducing the computa- 
tional cost of CFD analysis is to restrict atten- 
tion to linearized dynamics. The response of the 
linearized system about a nonlinear steady-state 
condition can be obtained via several methods. 
Some of these methods include the model or- 
der reduction of state-space models using various 
techniques 11 ’ 12 One method for building a lin- 
earized, low-order, frequency-domain model from 
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CFD analysis is to apply the exponential (Gaus- 
sian) pulse input . 13 This method is used to excite 
the aeroelastic system, one mode at a time, us- 
ing a broadband Gaussian pulse. The time-domain 
responses are transformed into the frequency do- 
main to obtain the frequency-domain generalized 
aerodynamic force (GAF) influence coefficient ma- 
trix. These linearized GAFs can then be used 
in standard linear aeroelastic analyses . 14 Raveh 
et al 15 applied this method while replacing the 
Gaussian pulse with step and impulse inputs. Re- 
cently, these time-domain, linearized impulse and 
step responses have been transformed directly into 
state-space form for use in other disciplines such 
as controls or optimization . 16 Guendal and Ces- 
nik 17 applied the Aerodynamic Impulse Response 
(AIR) technique, based on the Volterra theory, 
to the PM ARC aerodynamic panel code. The 
PM ARC/ AIR code was applied to a simplified 
High Altitude Long Endurance (HALE) aircraft 
for rapid linear and nonlinear aeroelastic analysis 
of the vehicle. Yet another method, different from 
the Volterra-based ROM approach, is the Proper 
Orthogonal Decomposition (POD) technique. The 
POD is a method that is used extensively at, sev- 
eral research organizations for the development 
of reduced-order models. A thorough review of 
POD research activities can be found in Beran and 
Silva . 18 A brief overview of the POD method and 
a representative result are included in this paper 
for completeness. In addition, a review of the is- 
sues involved in the development of reduced-order 
models for fluid-structure interaction problems is 
provided by Dowell and Hall . 19 A topic of re- 
cent interest is the potential development of hybrid 
POD/ Volterra methods. These hybrid techniques 
would combine the spatial resolution possible with 
POD methods with the low dimensionality and 
computat ional efficiency of Volterra methods. This 
is a very appealing concept that merits serious in- 
vestigation. 

A valuable and important characteristic of the 
Volterra theory of nonlinear systems is that the 
theory is well defined in the time and frequency do- 
mains for continuous- and discrete-time systems. 
In particular, this theory has found wide appli- 
cation in the field of nonlinear discrete-time sys- 
tems 20 and nonlinear digital filters for telecommu- 
nications and image processing. wl Howevei, ap- 
plication of nonlinear system theories, including 
Volterra theory, to modeling nonlinear unsteady 
aerodynamic responses has not been extensive. 
One approach for modeling unsteady transonic 
aerodynamic responses is Ueda and Dowell s~ ap- 
plication of describing functions, which is a har- 
monic balance technique involving one harmonic. 
Tobak and Pearson 23 apply the continuous-time 
Volterra concept of functionals to indicial (step) 


aerodynamic responses to compute nonlinear sta- 
bility derivatives. Jenkins 24 also investigates the 
determination of nonlinear aerodynamic indicial 
responses and nonlinear stability derivatives us- 
ing similar functional concepts, Stalford et aP 5 
develop Volterra models for simulating the behav- 
ior of a simplified nonlinear stall/post-stall aircraft 
model and the limit cycle oscillations of a simpli- 
fied wing-rock model. In particular, they establish 
a straightforward analytical procedure for deriving 
the Volterra kernels from known nonlinear func- 
tions. 

A particular response from a CFD code may 
provide information regarding the nonlinear aeroe- 
lastic behavior of a complex configuration due to 
a particular input at a particular flight condition. 

It does not, however, provide general information 
regarding the behavior of the configuration to a 
variation of the input, or the flight condition, or 
both. As a result, repeated use of the CFD code is 
necessary as input parameters and flight conditions 
are varied. A primary feature of the Volterra ap- 
proach is the ability to characterize a linearized or 
nonlinear system using a small number of CTD- 
code analyses. Once characterized (via step oi 
impulse responses of various orders), the functions 
can be implemented in a computationally-efficient 
convolution scheme for prediction of responses to 
arbitrary inputs without the costly repeated use of 
the CFD code of interest. Characterization of the 
nonlinear response to an arbitrary input via the 
Volterra theory requires identification of the non- 
linear Volterra kernels for a specified configuration 
and flight condition. Clearly, development and ap- 
plication of Volterra-based ROM concepts depend 
on the identification of the associated kernels for 
the problem of interest. 

The problem of Volterra kernel identification 
is addressed by many investigators^, including 
Rugh , 26 Clancy and Rugh, 2 ' Schetzen ,” 8 and more 
recently by Boyd, Tang, and Chua . 29 There are 
several ways of identifying Volterra kernels in the 
time and frequency domains that can be applied 
to continuous- or discrete-time systems. Tromp 
and Jenkins 30 use indicial (step) responses from 
a Navier-Stokes CFD code and a Laplace domain 
scheme to identify the first-order kernel of a pitch- 
oscillating airfoil. Rodriguez 31 generates realiza- 
tions of state-affine systems, which are related to 
discrete-time Volterra kernels, for aeroelastic anal- 
yses. Assuming high-frequency response, Silva 32 
introduces the concept of discrete-time, aerody- 
namic impulse responses, or kernels, for a rect- 
angular wing under linear (subsonic) and nonlin- 
ear (transonic) conditions. Silva 33 improves upon 
these results by extending the methodology to 
arbitrary input frequencies, resulting in the first 
identification of discrete-time impulse responses of 
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an aerodynamic system. It should be noted that 
owing to separation of the down wash input terms, 
Silva’s first approach had limited applicability for 
the identification of nonlinear Volterra kernels 33 a 
situation which has been resolved recently. 10 

In his dissertation, Silva 10 discusses the funda- 
mental differences between traditional, continuous- 
time theories and modern discrete-time formula- 
tions that allow the identification of discrete-time 
kernels. The discrete-time methods are then ap- 
plied to various nonlinear systems including a non- 
linear Riccati circuit, the viscous Burger’s equa- 
tion, an aeroelastic wing in transonic flow using a 
transonic small-disturbance code, and a supercriti- 
cal airfoil undergoing large plunge motions at tran- 
sonic conditions using a Navier-Stokes flow solver 
with the Spalart-AIlmaras turbulence model. Al- 
though the majority of the research mentioned 
thus far involving aeroelastic Volterra kernel iden- 
tification has been of a computational/analytical 
nature, the experimental identification of Volterra 
kernels for a nonlinear aeroelastic system has been 
performed. Kurdila et al 34 applied an efficient 
wavelet-based algorithm to the extraction of the 
nonlinear Volterra kernels of an aeroelastic system 
exhibiting limit cycle oscillations (LCO). There 
is increased interest in the development of these 
experimental techniques for use in various experi- 
mental settings. The identification of LCO during 
flight flutter tests is a case in point. 

One of the goals of the present paper is to pro- 
vide enough information to motivate readers to 
consult the references. The methods and results 
presented are not intended to be complete and all 
inclusive. Another goal of the present paperis to fa- 
miliarize the reader with NASA Langley Research 
Center's (LaRC s) vision for ROM research, unify- 
ing several of the methods and results presented. 
It should be stated that some of the research dis- 
cussed in the paper is funded by NASA while the 
other activities are purely of a cooperative nature. 
Development of ROMs using POD methods is cur- 
rently funded by the Air Force Research Labora- 
tory (AFRL) although a cooperative effort between 
the NASA LaRC The paper begins with a brief 
descript ion of the methods discussed: Volterra the- 
ory of nonlinear systems, state-space models from 
aerodynamic impulse responses, and POD. This is 
followed by presentation of various computational, 
analytical, and experimental results, including re- 
cent results based on the application of the POD 
method. The paper concludes with recommenda- 
tions for future research. 

Theoretical Background 

The research and development of reduced-order 
models (ROMs) for applications in nonlinear 
aeroservoelast-icity at the NASA Langley Research 



Fig. 1 Schematic representation of the 
Reduced-Order Modeling research program at 
the NASA Langley Research Center. 

Center can be categorized into three disciplines: 
computational fluid dynamics (CFD), system iden- 
tification, and control law design. A schematic 
representation of the components of this research 
program is presented as Figure 1. 

In addition to representing the vision for ROM 
research at the NASA LaRC, Figure 1 also serves 
as an organizational outline for this paper. The 
CFD portion consists of the development and ap- 
plication of techniques implemented into various 
CFD codes for the extraction of Volterra kernels. 
As will be discussed, the Volterra kernels represent 
linear and nonlinear functional ROMs for these 
CFD models. These functional ROMs can be used 
within a digital convolution setting to provide in- 
creased computational efficiency over the original 
CFD codes. These functional ROMs also can be 
transformed into the more traditional state-space 
models, suitable for use in modern control the- 
ory and optimization. The transformation of the 
Volterra kernels into state-space form is performed 
within the System Identification category. In ad- 
dition, the analytical derivation and the experi- 
mental identification of aeroelastic kernels will be 
discussed under this heading. The Control Law 
Design portion of Figure 1 is still in its early 
stages and is left as a topic for future publications. 
Research involving POD could easily fit within ei- 
ther the CFD or System Identification topics but 
is treated separately with various important and 
valuable references for description of details. 

Prior to discussion of the various results under 
the categories just described, a brief description of 
the various theoretical concepts follows. 

Volterra Theory 

We begin by reviewing key features of the 
Volterra theory, as applied to time-invariant, non- 
linear, continuous- and discrete-time systems. The 
literature on Volterra theory is significant, includ- 
ing several texts. 26 ’ 3:3-37 This section will con- 
centrate on time-domain Volterra formulations. 


4 OF 16 


International Forum on Aeroelasticity and Structural Dynamics Paper IFASD 2001-008 





consistent with the implied application to time- 
domain, computational aeroelasticity methods. 
Details regarding the foundations and applications 
of the frequency-domain Volterra theory can be 
found in several references. 26,29,36 39 Marzocca 
et al 38,39 present a thorough discussion of the 
frequency-domain Volterra theory with respect to 
nonlinear aeroelastic phenomena. This research is 
discussed in a subsequent section. 

We first consider time-invariant, nonlinear, 
continuous- time, systems followed by the applica- 
tion of Volterra theory to discrete-time systems 
(e.g., systems arising in CFD). Of interest is the re- 
sponse of the system about an initial state w(0) 
Wo due to an arbitrary input u{t) (we take u as a 
real, scalar input, such as pitch angle of an airfoil) 
for t > 0. As applied to these systems, Volterra 
theory yields the response 

w(/) = lio + / hi(< — t)u(t)cIt 
Jo 


+ f f h:>(< - r ht — T2)u(Ti)u(T'>)dridT2+ 
Jo Jo 



h n {t-Tu 


..,t-r n )u(ri)..u(r„)dn..dr„. 

( 1 ) 


The Volterra series in expression (1) contains three 
classes of terms. The first is the steady-state term 
satisfying the initial condition, lx 0 = W 0 . Next is 
the first response term, Jq lii(/ - t)u(t)(It , where 
hi is known as the first-order kernel (or the lin- 
ear/linearized unit impulse response). This term 
represents the convolution of the first-order kei- 
nel with the system input for times between 0 
and t. Lastly are the higher order terms involving 
the second-order kernel, I 12 , through the n th -order 
kernel, h n . The existence of these terms is an in- 
dication that the system is nonlinear. 40 

The convergence of the Volterra series is depen- 
dent on input magnitude and the degree of system 
nonlinearity. Boyd 41 shows that the convergence 
of the Volterra series cannot be guaranteed when 
the maximum value of the input exceeds a criti- 
cal value, which is system dependent. Of course, 
the issue of convergence is important, since the 
Volterra series must be truncated for analysis of 
practical systems. Silva 10,40 and Raveh et al. 
consider a weakly nonlinear formulation, where it 
is assumed that the Volterra series can be accu- 
rately truncated beyond the second-order term: 


w( / ) = ho + 



r)u(r)dr 


T 



h 2 (< 


n - T2)u(ri)u(rn)dTidT2. (2) 

5 


For linear systems, only the first-order kernel is 
non- trivial, and there are no limitations on input 
amplitude. 

Silva 10 derives the first- and second-order ker- 
nels, which are presented here in final form in 
terms of various response functions: 

h i{t) = 2wi(tl) - -w 2 (*l), (3) 

ho(<M2) = (wi(a,<2) - wi(tl) - w i(<2)) . 

( 4 ) 

In (3), wi {t 1) is the time response of the system to 
a unit impulse applied at time 0 and w 2 (21) is the 
time response of the system to an impulse of twice 
unit magnitude at time 0. If the system is linear, 
then w 2 = 2wo and hi = wq. If the system is 
nonlinear, then this ident ification of the fiist-order 
kernel captures an amplitude-dependent nonlinear 
effect. The identification of the second-order kernel 
is more demanding, since it is dependent on two 
parameters. Assuming t'2 > 1 1 in (4), wi(<2) is 
the response of the system to an impulse at time 
< 2 . 

Time is discretized with a set of time steps of 
equivalent size. Discrete time increments are in- 
dexed from 0 (time 0) to n (time t), and the 
evaluation of w at time n is denoted by w [??]. The 
convolution in discrete time is 
N 

w[n] = ho + hi[n — fe]w[fc] (fi) 

k = 0 

JV N 

- f h 2 [n — Aq, n — A* 2 ] u[Aq] u[Aro]. (6) 

k 1 zz0k 2 =0 

where N is the total time record of interest. 

It should be noted that an important conceptual 
breakthrough in the development and application 
of the discrete-time Volterra theory as a ROM 
technique is the distinction between a continuous- 
time unit impulse response and a discrete-time 
unit impulse response. 10,40 The continuous-time 
unit impulse response is an abstract function typi- 
cally defined as a function with an amplitude that 
reaches infinity while its width approaches zero but 
its integral is unity. This function is difficult, if 
not impossible, to apply in practical applications 
(i.e., discrete-time problems). The discrete-time 
unit impulse response (known as a unit sample 
response), on the other hand, is specifically de- 
signed for discrete-time (i.e., numerical) applica- 
tions. This function is defined as having a value 
of unity at one point in time and zero everywhere 
else. This is clearly a simpler function to imple- 
ment in a numerical setting. The proof of this and 
details regarding the very powerful unit sample re- 
sponse can be found in any modern text on digital 
signal processing. 42 
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The identification of linearized and nonlinear 
Volterra kernels is an essential step in the devel- 
opment of ROMs based on Volterra theory, but it 
is not the final step. Ultimately, these functional 
kernels can be transformed into linearized and 
nonlinear (bilinear) state-space systems that can 
be easily implemented into other disciplines such 
as controls and optimization. 10 26 Recently, lin- 
earized state-space models of an unsteady aerody- 
namic system have been developed while research 
into the development of nonlinear state-space mod- 
els continues. 16 

The frequency-domain version of the Volterra 
theory is simply the Fourier transform of the se- 
ries shown in (1). Therefore, the Fourier transform 
of the first-order kernel (for a linear system) is the 
frequency response function of the system. Higher- 
order kernels are Fourier transformed into higher- 
order frequency response functions, discussed in 
most of the references already mentioned. The 
primary benefit of these higher-order frequency 
response functions is that they provide informa- 
tion regarding the interaction of frequencies due 
to a nonlinear process. For example, bispectra 
(the frequency-domain version of the time-domain 
second-order kernel) have been used in the study of 
grid-generated turbulence to identify the nonlinear 
exchange of energy from one frequency to another. 
Linear concepts, by definition, cannot provide this 
type of information. In addition, some very in- 
teresting and fundamental applications using the 
frequency-domain Volterra theory 38,39 and exper- 
imental applications of Volterra methods 43,44 are 
providing new ‘‘windows" on the world of nonlin- 
ear aeroelasticity. 

Discrete-Time State-Space Models 

The basic formulation used in the development 
of state-space models of the unsteady aerodynamic 
system using the impulse responses (Volterra ker- 
nels) obtained directly from a CFD code is de- 
scribed in this subsection. The ability to generate 
state-space models of systems using impulse re- 
sponses was a primary motivation for the develop- 
ment of aerodynamic impulse response functions. 10 
Alt hough the method and results presented for this 
activity are linearized results, the long-term goal 
is the development of nonlinear state-space mod- 
els directly from the nonlinear Volterra kernels. 
The linearized results presented here are a starting 
point for the nonlinear state-space model develop- 
ment activity. Details can be found in the reference 
by Silva and Raveh. 16 

A finite-dimensional, discrete-time, linear, time- 
invariant dynamical system has the state-variable 
equations 

x[n + 1] = Ax[n] + Bti[n] (7) 


2/[n] = Cx[n] -f Du[n ] (8) 

where ,r is an m-dimensional state vector, u an p- 
dimensional control input, and y a r-dimensional 
output or measurement vector with n being the 
discrete time index. The transition matrix, A, 
characterizes the dynamics of the system. The 
goal of system realization is to generate constant 
matrices (A, B . C. and D) such that the output 
responses of a given system due to a particular set 
of inputs is reproduced by the discrete-time state- 
space system described above. 

For the system of (7) and (8), the time-domain 
values of the systems unit sample response are also 
known as the Markov parameters and are defined 
as 

Y[n] = CA n ~ 1 B (9) 

with B an m by p matrix and C a. q by m matrix. 
System realization techniques provide the constant 
matrices A , B . C and D using V(n). 

The Eigensystem Realization Algorithm (ERA) 
algorithm 45 46 begins by defining the generalized 
Hankel matrix consisting of the Markov parame- 
ters for all input/output combinations. The algo- 
rithm then uses the singular value decomposition 
(SVD) to compute the A , B , and C matrices. Of- 
ten. the direct feedthrough matrix, D , is nonzero 
whenever the initial values of the Markov parame- 
ters are nonzero. 

The ERA algorithm has been used successfully 
for the identification of several experimental struc- 
tural dynamic systems. Although the algorithm 
also has been used to extract damping and fre- 
quency information from CFD-generated aeroe- 
lastic transients (no published references), this 
research represents the first time that the ERA 
algorithm is applied to the development of un- 
steady aerodynamic state-space models using aero- 
dynamic pulse responses (Markov parameters). 
Additional details regarding the ERA algorithm 
and its numerous applications are discussed in the 
references provided. 

Proper Orthogonal Decomposition (POD) 

POD is a linear method for establishing an opti- 
mal basis, or modal decomposition, of an ensemble 
of continuous or discrete functions. The variables 
and parameters defined in this subsection are not 
related to variables from any previous subsection. 

Beran 18 provides an excellent summary of the 
POD method and its origins. Detailed deriva- 
tions of the POD and its properties are available 
elsewhere 4 ' 48 and not repeated herein. In this 
discussion of POD, M basis vectors are used to rep- 
resent deviations of w (t.) from a base solution. W 0 . 
These basis vectors are written as {eVe 2 , 
and are referred to by many names, including POD 
vectors, 49 empirical eigenfunctions 47 or, simply, 
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modes. 47,50 For the sake of brevity, we shall use 
the term “modes” to denote the POD basis vec- 
tors. Each mode is normalized such that e* e 1 = 1 
(i = 1, M), and computed in a manner to be 
described shortly. The modal decomposition of w 
using M modes, or w m, is given by 


AJ 


W m 


Wi^ = W 0 + *W, 


( 10 ) 


i = l 


where $ is an iV x M matrix containing the ordered 
set of modes, $ = [e 1 , e 2 , e M ] and w is an 
M- dimensional vector of modal amplitudes, w = 
[wi y W2, .... u?a/]. As a time-varying function, w is 
approximated by Wo + $w[t). 

As stated by Holmes et al. 47 “Linearity is the 
source of the [POD] method’s strengths as well as 
its limitations ...” The method is linear owing to 
the independence of the modes from the modal am- 
plitudes, thereby allowing for the straightforward 
construction of reduced-order equation sets from 
the full equation sets following mode computation. 

The POD modes are constructed by first com- 
puting samples, or snapshots, of system behavior 
(solutions at different instants in time for dynamic 
problems, or equilibrium solutions at different pa- 
rameter values for static problems) and storing 


hi 

force/plunge 



Fig. 2 Volterra kernels for CFD analysis 
of RAE airfoil:First-order kernel and effect of 
identification amplitudes. 



these samples in a snapshot matrix, S. This sam- 
pling process is often refered to as POD training. 
For now, we assume that M snapshots are collected 
and column- wise collocated into the N x M snap- 
shot matrix: S = [w 1 , w 2 , ..., w A/ ] . The basis 
provided by the POD, known as the Karhunen- 
Loeve 5152 (or I\-L) basis, has been shown to min- 
imize the error of the approximation of functions 

,1 1 r 1 • 47-49. 53, 54 

using these basis functions. 

In practice, fewer than M modes are retained 
to simulate system behavior. These are selected 
based on the size of the modal eigenvectors. Sim- 
ply put, the K-L basis for a subspace of dimension 
Mr < M is obtained by retaining the modes asso- 
ciated with the M r largest eigenvalues computed. 

The techniques described in Beran and Silva 18 
provide different means for obtaining reduced- 
order equation sets governing w(/) in the POD 
subspace. There are several methods for accom- 
plishing this including the Galerkin projection, 
“subspace” projection (for linear and nonlineal 
systems), and collocation. These methods are dis- 
cussed in detail by Beran and Silva and will not 
be repeated here. In addition, recent results by 
Kim and Bussoletti 55 and Kim 50 are indicative of 
the efficiency and suitability of these techniques. 
However, due to the importance and value of POD 
techniques, a sample result is provided in a sub- 
sequent section in order to familiarize the reader 
with the appropriate references concerning POD 
research. 


Fig. 3 Volterra kernels for CFD analysis 
of RAE airfoil: First five components of the 
second-order kernel for plunge amplitude of 0.1. 

Results 

The results discussed in this section address the 
top two disciplines (circles) of Figure 1. The hist 
subsection discusses the approximation of CFD 
results using ROMs and the second subsection dis- 
cusses results from system identification studies. 

Approximating CFD Results with ROMs 

This section describes some results obtained 
in the identification of CFD-based (time-domain) 
Volterra kernels. The first result presented is from 
Silva 10 for a supercritical airfoil at a transonic 
Mach number. 

First- and second-order kernels for the Navier- 
Stokes solution (with Spalart-Allmaras tuibulence 
model) of an RAE airfoil undergoing plunging 
motions at a transonic Mach number using the 
CFL3D code 57 are presented in Figures 2 and 
3. Shown in Figure 2 are two sets of first- 
order kernels due to two different sets of excita- 
tion amplitudes. Recall that the first-order kernel 
is computed using (3) and is the result of two 
pulse responses, one at a particular amplitude and 
the second at double the first amplitude. One 
of the first-order kernels shown in Figure 2 was 
computed using excitation plunge amplitudes of 
w - 0.01 and w — 0.02, where u? is a fraction 
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of the chord of the airfoil. The other first-order 
kernel was computed using excitation plunge am- 
plitudes of w = 0.1 and w = 0.2. It is clear that the 
two kernels are not linearly related, demonstrating 
how the first-order kernel can capture amplitude- 
dependent nonlinear effects. 

Shown in Figure 3 are five components of the 
second-order kernel for this case. The second-order 
kernel is more complicated because it is a two- 
dimensional function of time. The second-order 
kernel is presented as a family of functions in Fig- 
ure 3 for simplicity. The important point to 
be made is that this kernel is readily generated 
and its relatively smaller values (compared with 
the first-order kernel) and its rapid convergence 
indicate a small (but not negligible) level of non- 
linearity present. This information may be used 
to determine if t he first-order kernel is sufficient to 
capture the dominant nonlinear effects. This point 
is demonstrated in Figure 4. 

Figure 4 is a comparison of plunge responses 
for three different plunge amplitudes for the same 
configuration. Specifically, three comparisons are 
made between the full CFD solution due to a si- 
nusoidal plunging motion (labelled in the figure 
as “actual nonlinear'’) and that obtained using 
the first-order kernel from Figure 2 due to the 
larger excitation amplitudes. As can be seen, the 
plunge response obtained using the Volterra first- 
order kernel compares almost identically with the 
response obtained from the full CFD solution for 
the two smaller amplitude responses. In fact, the 
two curves for these responses, are almost not dis- 
cern able in Figure 4. The comparison for the 
largest amplitude response is very good as well, 
with a slight but noticeable difference between the 
two results. The nonlinearity of the large- ampli- 
tude plunge responses is confirmed by linearly scal- 
ing the smallest amplitude (i.e., linear) response 
which, as shown in Figure 4, cannot capture the 
amplitude-dependent nonlinear plunge dynamics 
seen at the larger amplitude. 

The turnaround time (“wallclock”) for the full 
CFD solution was on the order of a day whereas 
the Volterra first-order solution was computed on 
a workstation in 30 seconds using digital convolu- 
tion. The initial cost of computing the first-order 
kernel was trivial as well due to the rapid conver- 
gence of the pulse responses. In fact, since each 
pulse (unit and double amplitudes) goes to zero in 
less than 100 time steps, the responses were gener- 
ated using a debugging mode option available on 
the supercomputer system used. Using this option, 
computations requiring less than 300 time steps 
are executed immediately, intended for debugging 
purposes. As a result, each pulse was computed 
within five minutes, resulting in a first -order ker- 
nel that was computed in about ten minutes. Of 



Fig. 4 Comparison of actual nonlinear and 
first-order Volterra responses for three differ- 
ent plunge motions and a linear response for 
the largest motion. 

course, once t lie kernel has been computed, it can 
be used to predict the response to any arbitrary 
input (steady, any and all frequencies, random) of 
arbitrary length via digital convolution on a work- 
station. Using this method, there is no need for the 
repeated, and costly, execution of the CFD code for 
different inputs. 

Raveh, Levy and Karpel have recently applied 
the Volterra- based ROM approach to analysis 
of the AGARD 445.6 wing. 15 They simulate 
the flow field around the wing using the EZNSS 
Euier/Navier-Stokes code. 58 This code provides a 
choice between two implicit algorithms, the Beam- 
Warming algorithm 59 and the partially flux-vector 
splitting algorithm of Steger et al. 60 Grid genera- 
tion and inter-grid connectivity are handled using 
the Chimera approach. 61 The code was enhanced 
with an elastic capability to compute trimmed ma- 
neuvers of elastic aircraft. 58 For the CFD compu- 
tations, the flow field around the wing was evalu- 
ated on a C-H type grid, with 193 points in the 
chord wise direction along the wing and its wake, 
65 grid points in the spanwise direction, and 41 
grid points along the normal direction. 

A process of mode- by-mode excitation, dis- 
cussed previously, was performed for this wing 
using four elastic modes at a Mach number of 
0.96. The mode-by-mode excitation technique pro- 
vides the unsteady time-domain generalized aero- 
dynamic forces (GAFs) in all four modes due to 
an excitation of one of the modes. In this fash- 
ion, a matrix of four- by-four functions (sixteen 
total) is developed. Two sets of excitation inputs 
were used: the discrete-time pulse input and the 
discrete-time step input. The cost of computing 
these functions is minimal due to the rapid con- 
vergence of these functions and it consists of only 
four code executions. Once these functions were 
defined, several full CFD solutions were generated 
that were due to various sinusoidal inputs at vari- 
ous frequencies. Shown in Figure 5 is just one of 
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Fig. 5 Comparison of actual and first-order 
Vo l terra responses to sinusoidal excitation at 5 
Hz (Mach 0.96). 

these results for a 5 Hz input frequency, comparing 
the result obtained from the full CFD solution to 
that obtained via convolution of the step or pulse 
responses with a 5 Hz sinusoid. As can be seen, 
the comparison is exact to plotting accuracy for 
most of the responses. The full CFD solution, con- 
sisting of 8000 iterations required approximately 
24 hours on an SGI Origin 2000 computer with 4 
CPUs. By comparison, the Vol terra- based ROM 
response shown required about a minute. Even 
including the upfront cost of computing the pulse 
(or step responses), the computational cost savings 
are significant. More importantly, the same pulse 
(or step) functions can now be used to predict the 
response of the aeroelastic system to any aibitiaij 
input of any length. 

As a validation, a comparison of this approach 
to another result available in the liteiatuie is pie- 
sented as Figure 7. Shown in Figure 7 is a com- 
parison of linear and nonlinear GAFs for the first 
two modes of the AGARD 445.6 Aeroelastic \\ ing 
at a Mach number of 0.96. Nonlinear GAFs refers 
to the GAFs computed using the Volterra pulse- 
response technique about a nonlinear steady-state 
value by exciting one mode at a time and obtaining 
the resultant responses in the other modes. The 
CFD results are compared with those using the 
ZAERO code for a purely linear case. Frequency- 
domain values were obtained by performing a con- 
volution of several frequencies of interest with the 
computed CFD-based pulse responses. As can be 
seen, the comparison is reasonable and shows the 
small (but not negligible) nonlinear aerodynamic 
effects identified using the Volterra pulse-response 
technique. 

Additional computational applications of the 
Vol terra- based ROM technique include the appli- 
cation to a High Altitude Long Endurance (HALE) 
wing using a panel method for the aerodynamics. 
The simplified wing is referred to as the Simple 
High Aspect Ratio (SHAR) Wing. The particu- 
lar panel method selected for this research was the 
PM ARC aerodynamic panel code, a code devel- 
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Fig. 6 Pitch responses of the SHAR Wing to 
a prescribed pitch command. 

oped in the late 1980s at NASA Ames Research 
Center. 62 This code is a low-order panel method, 
meaning that the source and doublet singulari- 
ties are considered to be constant on each panel. 
PM ARC is one of the most advanced panel meth- 
ods available, featuring an advanced time-stepping 
relaxed wake, internal flow capabilities, jet plume 
modeling, and rudimentary unsteady capabilities. 

The implementation of the Volterra-based ROM 
technique into the PM ARC code is referred to 
as the Aerodynamic Impulse Response (AIR) 
technique. A particularly useful feature of 
PM ARC/ AIR is the ability to predict the response 
to an arbitrary input. For the purpose of testing 
this feature, the following Fourier sine series with 
three components was used: 

9[t] = 3° sin(27rt) + 3° sin(4^) + 3° sin(6^<) (11) 

The quantity theta [t] is a signal with three fre- 
quencies: 1 Hz, 2 Hz, and 3 Hz. The pitch response 
of the SHAR wing is shown in Figure 6. The re- 
sponse was calculated with the AIR method using 
impulse responses generated with three time steps, 
At = 0.01 seconds, 0.015 seconds, and 0.02 sec- 
onds, and with the full solution from the PMARC 
code at At = 0.015 seconds. All three time step 
choices produce a fairly accurate prediction of the 
response. These results indicate that the method 
is a useful one for predicting the response to arbi- 
trary inputs, especially those that may not be as 
easy to implement within the PMARC code. This 
is a significant improvement over the capabilities 
of the original PMARC code. 

System Identification 

In this section, some of the results classified un- 
der the svstem identification portion of Figure 1 
are discussed. These results include the tranfor- 
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mat ion of Volterra kernels into state-space form, 
and analytical and experimental identification of 
Volterra kernels for various systems. The state- 
space resuts follow the method described in an 
earlier subsection. The analytical identification 
of Volterra kernels consists of two parts: 1) the 
derivation of compressible aerodynamic indicia! 
functions for subsonic, supersonic, and hypersonic 
Mach number regimes and 2) the derivation of 
frequency-domain Volterra kernels for aeroelastic 
systems with structural nonlinearities. It should 
be pointed out that the relationship between in- 
dicia! functions and Volterra kernels is the same 
as the relationship between a step and impulse re- 
sponses: one is the derivative of the other. The 
experimental identification of Volterra kernels con- 
sists of the experimental extraction of second-order 
Volterra kernels for an aeroelastic system under- 
going limit cycle oscillations (LCO). Results from 
these categories are now discussed. 

In an earlier subsection, the generation of time- 
domain GAFs using a CFD code was presented. 
Validation of these functions was performed in the 
time and frequency domains. Time-domain vali- 
dation consisted of a comparison between the full 
CFD solution and the solution obtained via con- 
volution for forced harmonic responses at several 
frequencies. Frequency-domain validation was per- 
formed by Fourier transforming the time-domain 
GAFs to the frequency domain and comparing 
these with results available in the literature for 
that configuration. Because traditional flutter 
analyses are performed in the frequency domain 
using frequency domain GAFs, the standard ap- 
proach for using linearized GAFs from a CFD 
code is to transform the time- domain GAFs into 
the frequency domain and use the standard flutter 
analysis routines. Then, if a time-domain model of 
the GAFs is needed for ASE analyses (i.e., state- 
space form), then rational function approximations 
(RFAs) are employed to transform the frequency- 
domain GAFs back into the time domain. But 
rather than transforming the time-domain GAFs 
into the frequency domain only to transform them 
back into the time domain, discrete-time, state- 
space systems can be created using the Volterra 
kernels (time-domain GAFs) directly. 16 

Using the ERA method previously discussed, a 
32nd-order state-space model of the AGARD 445.6 
Aeroelastic Wing was generated. This state-space 
model consists of four inputs and four outputs, one 
for each of the four structural modes. Then, in 
order to validate the accuracy of the state-space 
model, the frequency content of the state-space 
model is compared to the frequency content of 
the Volterra kernels extracted from the CFD code. 
Presented in Figure 8 is a comparison of the fre- 
quency response for the CFD solution (kernels) 



Fig. 7 Comparison of linear and nonlinear 
GAFs for the first two wing inodes. 



Fig. 8 Comparison of frequency responses for 
AGARD 445.6 Aeroelastic Wing. CFD (Solid 
Line), State-Space (Dotted Line). 


versus the frequency response for the state-space 
system for the AGARD 445.6 Aeroelastic Wing. 
Presented in the figure are the responses (output) 
for all four modes due to an input in the first mode. 
The frequency responses of the pulses computed 
directly from the CFD code (solid lines) compare 
very well with the frequency responses obtained 
from 32nd-order state-space system generated to 
model this system. The responses due to inputs 
in the second, third, and fourth modes are just as 
good as those shown in Figure 8, but are not pre- 
sented here for brevity. The important point to be 
made is that time-domain kernels, extracted from 
a CFD code, can now be transformed directly into 
state-space form, amenable for use in modern con- 
trol theory and optimization. 

In terms of analytically-derived kernels (indi- 
cia! functions), Marzocca et al 6 have developed 
a unified approach, based upon the use of aero- 
dynamic indicia! functions, that enables the de- 
termination of the aerodynamic lift and moment 
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Fig. 9 Comparison of aeroelastic responses due 
to an external blast load. 


in subsonic, supersonic, and hypersonic compress- 
ible flight speed regimes for three-dimensional 
wings. These indicia] functions are the compress- 
ible, three-dimensional extensions of the classical 
incompressible, two-dimensional indicial functions, 
such as Wagner's function. As a result, these in- 
dicial aerodynamic functions can now be used to 
perform rapid evaluations of preliminary designs 
over a wide range of Mach number regimes and for 
three-dimensional wings. 

A sample result from this research for a swept 
wing in the supersonic and hypersonic flight speed 
regimes is presented as Figure 9. Shown in Figure 
9 is a comparison of the aeroelastic responses due 
to a blast load using the method of Marzocca et 
al 6 and piston theory for various Mach numbers. 
As can be seen, the comparisons are excellent, vali- 
dating this methodology for this speed range. Sim- 
ilar comparisons for subsonic speed regimes can 
be found in the stated reference. A fundamen- 
tal contribution of this research is that it defines 
the relationship between aerodynamic indicial re- 
sponses and aerodynamic Volterra kernels. This 
functional relationship can then be used to deiive 
nonlinear aeroelastic Volterra kernels. 

Analytical derivation of frequency-domain, non- 
linear aeroelastic Volterra kernels is presented by 
Marzocca et al. 38,39 Aeroelastic systems with one 
(plunge) and two (plunge and pitch) degrees of 
freedom with linear aerodynamic equations but 
with nonlinear structural parameters were defined. 
The nonlinear structural parameters consisted of 
quadratic and cubic stiffness and damping terms. 
Analytical derivation of the Higher-Order Fie- 
quency Response Functions (HO-FRF), which is 
another term for frequency-domain Volterra ker- 
nels, was performed on these nonlinear aeroelastic 
systems. Evaluation of the HO-FRFs was per- 
formed by comparing the responses obtained using 
these HO-FRFs with the responses obtained via di- 
rect simulation of the systems. 

Presented in Figure 10 are the results of a con- 
vergence study involving the first three kernels for 


First Approximation 

Second Approximation 

Third Approximation 

— Exact Response (Num. Int.) 




Fig. 10 Convergence study involving the first 
three kernels and comparison with the exact 
nonlinear aeroelastic response to a 1). 1-cosine 
gust and a 2). triangular external blast load, as 
shown in insets. 

the one-degree-of-freedom system and comparison 
with the u exact” nonlinear aerolastic response to 
a (a) 1-cosine gust pulse and a (b) triangular blast 
load as shown in the insets of the figure. These 
time-domain results were obtained by computing 
the inverse Fourier transform of the derived HO- 
FRFs. This demonstrates the rapid and accurate 
convergence of the Volterra kernel approach for 
these types of problems. 

For the two-degree-of- freedom system, a three- 
dimensional plot of the second-order HO-FRF 
is presented as Figure 11, where the theoret- 
ical symmetry about the diagonal elements of 
these Volterra kernels is verified. The importance 
of these results is that these analytically-derived 
functions provide a closed-form nonlinear solution 
of simple nonlinear aeroelastic systems which can 
be used for the validation of computational tech- 
niques. In addition, enhanced understanding of 
these HO-FRFs can provide significant insight into 
the nonlinear dynamics of a system that would not 
be visible via standard linear (or linearized) pro- 
cesses. 

A discussion of the experimental results, as per- 
formed by Kurdila et al 43 and Prazenica et al 44 
now follows. First of all, it is important to mention 
that, historically, nonlinear Volterra series have 
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Fig. 11 Three-dimensional and contour plots 
of second-order aeroelastic kernel. 


not seen widespread use in system synthesis pre- 
cisely because of the high dimensionality of the 
higher order, nonlinear terms. This is true from 
experimental, computational, and analytical per- 
spectives. However, recent work by researchers in 
multiresolution analysis of the Volterra kernels has 
shown that the dimensionality of the higher order 
terms can be significantly reduced. This reduction 
is due to the fact that wavelet and multiresolu- 
tion analysis have shown considerable promise for 
the compression of signals, images, and, in partic- 
ular, some integral operators. Numerical studies 
of the theoretical derivations are carried out using 
experimental pitch and plunge response data from 
the Texas A & M University’s (TAMU) Nonlinear 
Aeroelastic Testbed (NAT). 

The TAMU-NAT, as used for this study, was 
comprised of a NACA0012 airfoil with a flap 
mounted on a specially-designed carriage capable 
of undergoing large amplitude pitch and plunge 
motions. Aside from the potential nonlinear aero- 
dynamic effects int roduced by the stall of the airfoil 
at large angles of attack, this experimental setup 
results in an aeroelastic system with, primarily, 
nonlinear structural effects. The nonlinearity in 
the system is introduced by the appropriate vari- 
ation of the various parameters of the pitch-and- 
plunge carriage such as cams and springs. 

A wavelet-based kernel identification algorithm 
was developed to extract Volterra kernels from 
experimental data. This algorithm was then ap- 
plied to the experiment involving the TAMU-NAT. 
The information used by the wavelet-based ker- 
nel identification algorithm consisted of the flap 
deflection (input) and pitch angle (output) as func- 
tions of time. Variation of the flow velocity of 
the TAMU-NAT ranged between 23 and 24 m/s, 



Time, sec 

Fig. 12 Comparison of predicted output from 
model identification and experimental output. 


which classifies this system as a time- varying or 
noil-stationary system. Strictly speaking, this clas- 
sification is not consistent with the theoretical de- 
velopments for stationary nonlinear systems. How- 
ever, this flow velocity variation was seen to be 
small enough to classify the system as “nearly sta- 
tionary”. In addition, windowed sampling of the 
data was performed in order to further reduce the 
non-stationarity of the system. 

Presented in Figure 12 is a comparison of the 
predicted output from the model identification and 
the experimental output for the airfoil exhibiting 
limit cycle oscillations (LCO). As can be seen, 
the two results are indistinguishable, verifying the 
accuracy of the identified model. In addition, Fig- 
ures 13 and 14 are second-order kernels identified 
for this system at different, samplings (temporal 
windows) during the experiment. A noticeable 
difference exists between the second-order kernel 
of Figure 13 and that of Figure 14, which is 
indicative of a modification of the nonlinear sys- 
tem. Clearly, some event took place during the 
test that altered, in some sense, the nonlinear sys- 
tem. This change is also visible in Figure 12 as 
a noticeable change with time in the nature of the 
time series of the output. This change in the non- 
linear dynamics of the system is not due to flow 
velocity variations but the exact cause is not yet 
known. These results demonstrate and validate 
two fundamental properties of Volterra kernels: 1) 
the existence of the higher-order kernels is an in- 
dication of a nonlinear process or system, and 2) 
alterations of the nonlinear dynamics of a system 
are reflected as alterations of these kernels. Finally, 
the identification of the second-order kernels was 
performed using as few as 14 wavelets, significantly 
reducing the computational burden as compared 
to more traditional Volterra kernel identification 
techniques. 
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Fig. 13 Vol terra second-order kernel for initial 
portion of output data. 
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Fig. 14 Vol terra second-order kernel for latter 
portion of output data. 


POD Results 

To assess the applicability of POD-based ROMs 
to differential equations that exhibit limit-cycle os- 
cillation (LCO), Beran et al. 63 computed solutions 
with the subspace projection technique of a tubu- 
lar reactor, known to experience LCO" 40 The 
governing equations are 

— = Lwi - wiQ(wi), (12) 

dt 


— = Lwn - j3i (u-2 - &) + fomQlu’i), (13) 

di 


_ i d 2 
L = Ye 8x- 


where Pe, /?i, {3 2 , Pj f* aie specified pa- 
rameters. Equations (12) and (13) describe con- 
vection, diffusion and reaction (CDR) within the 
reactor, and are referred to as the CDR equations. 


The variables uq and w 2 represent concentration 
and temperature, respectively, and the parameter 
(the Damkohler number) determines the abil- 
ity of the CDR equations to sustain LCO. The 
spatial domain is normalized; boundary conditions 
are applied at r = 0 and x = 1* Following 
spatial discretization of the equations and spec- 
ification of suitable initial conditions, which are 
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Fig. 15 Comparison of full-system and ROM 
solutions for an 8-inode ROM trained at n — 
0.16. Equilibrium and LCO branches. 


described elsewhere, 63 the CDR equations take the 
form ^ = R(w;p). 

The CDR system experiences a supercritical 
Hopf bifurcation at f = n = 0.1 6504, 63 which 
is accurately predicted by a POD-based ROM. 
The stability properties of the CDR system are 
shown in 15, where it is seen that stability of 
the equilibrium branch is lost beyond the bifur- 
cation point. Solutions are characterized by the 
maximum value of temperature computed over the 
domain, T max . The ROM is developed by sampling 
the CDR system as it evolves towards steady-state 
(0 < t < 2.5) for a value of /.< leading to system 
stability and starting with //° = 0.16 and P° < ^ • 
Following the procedure described above, 8 modes 
are computed and retained, representing a 15-fold 
reduction in problem size. Solutions of the full sys- 
tem are explicitly computed via time integration. 

Equilibrium solutions of the full system and the 
ROM are observed to be in excellent agreement. 
As seen in Figure 15, agreement is nearly exact 
— 0.16, where the POD is constructed, and 
is excellent for the remaining values of ji shown. 
Beyond the Hopf point, LCO amplitude is well pre- 
dieted with the ROM. 

The critical value of (x at which the CDR sys- 
tem loses stability is also very accurately predicted 
using reduced-order modeling. Stability loss is 
observed at [x = 0.16503, nearly the same value 
predicted with the full-system equations (/i = 
0.16504). 

It is clear that POD and related methods pro- 
vide insight into the dominant dynamics of a given 
system. Future research may focus on the coupling 
of Vol terra methods with POD methods for a hy- 
brid approach. 


13 OF 16 


International 


Forum on Aeroelasticity and Structural Dynamics Paper IFASD 2001-008 




Concluding Remarks 

The development of ROMs based on the time- 
domain and frequency-domain versions of the 
Vol terra theory of nonlinear systems has been de- 
scribed, including continuous- and discrete-time 
versions of the theory. The basic objective of the 
t heory is the identification of linearized and nonlin- 
ear kernel functions that capt ure the dominant re- 
sponse features of a nonlinear system. The method 
provides a very natural and intuitive extension of 
well-understood linear phenomena into the nonlin- 
ear domain. The results presented highlight two 
primary advantages of the theory: improved com- 
putational efficiency and insight into the nonlinear 
dynamics of the system of interest. 

The status of the Vo It err a- based ROM approach 
can be summarized as follows. The method has 
been used to show that discrete-time concepts, in- 
deed digital signal processing concepts such as unit 
pulses and step inputs, are directly applicable to 
CFD codes. The method has also been shown to 
be a higher-level generalization of the st andard lin- 
ear methods in use today. In addition, the nature 
of the method is such that it requires minimal 
modification to the CFD code of interest. Most 
unsteady aerodynamic or aeroelastic CFD codes 
already have various excitation inputs (e.g., sinu- 
soidal) and extension to a Volterra-based ROM 
approach simply involves adding a pulse (or step) 
input to the suite of available inputs - the CFD 
code itself remains unchanged. 

As for POD- based ROM developments (and re- 
lated methods), it is clear by the quantity and 
quality of the research effort devoted to this ap- 
proach that a great deal of insight can be gained 
into the system being analyzed. For example, 
a POD-based ROM of a discretized convection- 
diffusion-reaction (CDR) system was described 
and shown capable of determining a variety of 
important characteristics of the nonlinear system, 
including nonlinear static behavior, bifurcation to 
limit-cycle behavior, and sensitivity to changes in 
system parameters. The CDR problem serves as 
an analog for the study of the aeroelastic proper- 
ties of a wing, including static analysis, dynamic 
analysis, bifurcation analysis, and sensitivity anal- 
ysis. Several of the POD references mentioned thus 
far provide detailed discussions regarding future 
extensions and applications of the method. 

Future Directions 

As for the challenges associated with the 
Volterra-based ROM approach, there is much work 
to be done. An important issue that needs to 
be addressed is the issue of modal superposition 
with respect to nonlinear effects. Although it is 
clear that a mode- by- mode excitation is a lin- 
earization of the aeroelastic process, it is important 


to understand the limitations of this approach. 
In addition, work continues on the development 
of a technique that provides simultaneous exci- 
tation to all modes, eliminating the linearization 
issue. Linearized state-space models have been 
developed using the CFD-based pulse responses 
that can be incorporated directly into control sys- 
tem analysis, for example. These state-space ma- 
trices also sidestep the need to transform time- 
domain CFD loads into the frequency-domain only 
to transform the frequency-domain loads back into 
the time domain via rational function approxima- 
tions. Using the Volterra approach, time-domain 
CFD-based information goes directly into creat- 
ing time-domain state-space matrices, a more ef- 
ficient process. But the ultimate challenge lies 
in the creation of nonlinear (bilinear) state-space 
matrices which are mathematically related to the 
higher-order Volterra kernels. Some work has been 
done in this area, but there is significantly more 
work that needs to be done. With respect to the 
frequency-domain version of the theory, additional 
research is needed to develop precise interpreta- 
tions of the higher-order spectra as they apply 
to nonlinear aeroelastic systems. Additional re- 
search in the analytical derivation of these func- 
tions as well as experimental application of these 
techniques promises to provide great insights to all 
these questions. 
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